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Motivated by recent proposals on strain-engineering of graphene electronic circuits we calculate 
conductivity, shot-noise and the density of states in periodically deformed graphene. We provide the 
solution to the Dirac-Kronig-Penney model, which describes the phase-coherent transport in clean 
monolayer samples with an one-dimensional modulation of the strain and the electrostatic potentials. 
We compare the exact results to a qualitative band-structure analysis. We find that periodic strains 
induce large pseudo-gaps and suppress charge transport in the direction of strain modulation. The 
strain-induced minima in the gate-voltage dependence of the conductivity characterize the quality 
of graphene superstructures. The effect is especiaUy strong if the variation of inter-atomic distance 
exceeds the value /I, where a is the lattice spacing of free graphene and £ is the period of the 
superlattice. A similar effect induced by a periodic electrostatic potential is weakened due to Klein 
tunnelling. 



PACS numbers: 73.23.-b, 73.22.Pr, 73.21. Cd 

I. INTRODUCTION 

Graphene is recognized as the only two-dimensional 
crystal that withstands large deformations and high 
temperatures and is readily integrated into electronic 
circuits ji"— The high electron mobility and unique spec- 
tral characteristics single out graphene as a promising key 
element for electronic and optoelectronic applications 

The well-known Dirac dispersion of electronic ex- 
citations in graphene is responsible for a num- 
ber of interesting analogies with relativistic quantum 
electrodynamicSfiS but prevents the use of graphene in 
field-effect transistors due to the absence of energy band- 
gaps. ^iiiS The strain-engineering of graphene circuits 
has been recently suggested as a means to bypass this 
difficultyJ^ii^ In the Dirac-fermion picture of graphene, 
the strain gradient is equivalent to the presence of a 
pseudo-magnetic field, which can be manipulated to in- 
duce a zero-field quantum Hall effect and a topological 
isolator state.J^ The variation of strain may also induce 
substantial variation in the on-site electron energies, how- 
ever, this effect is reduced by screeningJ^ 

The experimental realizations of strained graphene 
confirm its surprisingly high elasticity and pave the 
way for the development of strain-engineered graphene 
electronics j^iii Different techniques have been proposed 
to produce graphene samples with controlled periodical 
variations of straiuiJ^^-Stt First transport measurements of 
graphene superstructures have been already reported.^^ 
Below we analyze the phase-coherent charge transport 
in periodically strained graphene samples and propose a 
way to characterize the quality of graphene superstruc- 
tures on the basis of their transport properties. 

More specifically we use the exact solution of the Dirac- 
Kronig-Penney model to calculate density of states, con- 
ductance, and shot noise in transport through finite size 
graphene samples with periodic potentials. The scatter- 




FIG. 1: (Color onhne) Schematic illustration of the graphene 
setup with metal leads for a; < and x > L. The angle 6 
specifies the orientation of the honeycomb lattice with respect 
to the transport direction x. 



ing off the metal leads is taken into account in all quan- 
tities. The position, the width, and the shape of the con- 
ductance minimum as well as the shot noise maximum 
associated with the periodic superstructure are analyzed 
in detail. 

This paper is organized as follows. In Section |TT] the 
effective Dirac Hamiltonian for deformed graphene is de- 
rived. The scattering approach to transport is briefly de- 
scribed in Section Hm Section ITVl is devoted to the exact 
solution of the Dirac-Kronig-Penney model for graphene 
with one-dimensional modulations of strain in transport 
direction. The generalization of this model is discussed 
at the end of the Section IIVI The results for conduc- 
tance, shot-noise, and the density of states are quali- 
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tatively understood on the basis of the band structure 
analysis, which is presented in Appendix |^ In Section 
lIVDI we discuss transport in the direction perpendicular 
to the strain modulation. We summarize our results in 
the concluding Section IVl 



II. EFFECTIVE HAMILTONIAN 

The electromechanical coupling in deformed graphene 
membranes has been investigated theoretically by many 
authors^Sr^^ following earlier publications on carbon 
nanotubes.— The deformations affect the hopping in- 
tegrals in the tight-binding description of graphene in two 
distinct ways: by changing the distance between carbon 
atoms and by tilting the electronic p^-orbitals, which are 
responsible for conduction. However, for most cases of 
interest, the in-plane strains play the major role in de- 
termining the electronic properties of deformed graphene, 
while the tilting can be neglected. 

The effects of in-plane strain are well captured by the 
tight-binding hamiltonian 
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(1) 



R a = l 



where the summation runs over the atomic positions, R, 
of a honeycomb lattice, the annihilation operators ur 
and bR refer to the occupation of Pz-orbitals at the two 
non-equivalent positions, A and B, of the unit cell, and 
the three vectors Sa shown in Fig. [T]are directed from a 
B-atom to its three nearest neighbors. 

If the strain varies smoothly on atomic distances, the 
deviation of the hopping integral from its unperturbed 
value t « 2.7 eV in a perfect crystal can be parameterized 
as 



St^iR) = {l3t/a^)dluiR)5o., 



(2) 



where f3 = —d\nt/d\na w 2 and u{R) is the strain ten- 
sor of the graphene membrane. 

The dimensionless strain tensor describes a change in 
a metric that can be expressed as 



(3) 



where the summation over the spatial indices i,j — x,y 
is assumed. The length elements c?£q = Sijds^ds'-' and 
dp = gijds''ds^ correspond to the metric Sij in the flat 
space and to the local metric gij in the membrane, re- 
spectively. In subsequent formulas we let conventionally 
s'^ = X and — y. The graphene crystal withstands 
very large internal strains so that the values of the strain 
tensor elements may reach 20% 42 

Even though the strain tensor makes no reference to 
the crystal structure of graphene, the lattice symmetry 
is entering the hopping integral in Eq. ^ due to the 
vectors (5q, a = 1,2,3. The slow variation of u{R) on 
the scale of the lattice spacing justifies the effective mass 



approximation, which is formulated in a generic reference 
frame shown in Fig. [1] In this frame wc find 



-cos(0 + 0„) 



1 



27ra/3, (4) 



while the positions of the two non-equivalent Dirac points 
in the reciprocal space are given by 



„ 47r /cos6l 



(5) 



Using the Fourier ansatz 

2 2 

aR = Y,e^'''''a,{R), 6r = ^ e*-f^=%(H), (6) 

s=l s=l 

we obtain, to the leading order in spatial gradients, the 
effective model 

H = j d^R [-ifivaV + T^crAg^R)] ^ r, (7) 

where hv — 3to/2, is the Pauli matrix in the valley 
space and the operators are arranged into the four-spinor 
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(8) 



The vector field Ag = {Agx, Agy) is real and its compo- 
nents are related to the strain tensor in a simple way. 



Agx - iAgy = {hvl3/a)e''\u^.j: + 2m, 



xy 



)■ (9) 



In the derivation of the eff'ective Hamiltonian ([7]) we ne- 
glected the velocity renormalization, which would appear 
as a correction of the order Sta to the prefactor of the 
spatial gradients.?- 

Unlike the usual vector potential, the vector field Ag 
preserves the time reversal symmetry of the Hamiltonian 
(l7|). It also exhibits the discrete rotational invariance 
of the honeycomb lattice. Indeed, the Hamiltonian ([7]) 
remains invariant under the rotation through the angle 
= 27r/3, while the rotation through the angle 9 = n/i 
is equivalent to an interchange of valleys. 

The eff'ect of a constant uniaxial strain has been stud- 
ied both theoreticallj^i23. and experimentally^i^i and 
falls beyond the scope of our consideration. If the con- 
stant uniaxial strain (in x direction) exceeds a certain 
critical value, the Dirac points merge and a band gap 
opens. It is predicted"^^ that, for a crystal expanded uni- 
formly in the zigzag direction (0 = 0), the critical expan- 
sion takes on its minimal value (~ 23%), while uniaxial 
strain in armchair direction, 9 = 7r/2, never generates a 
gap. These predictions await experimental verification. 

The strain also induces an electrostatic potential due 
to a change of the on-site energies of the tight-binding 
model. This effect leads to the appearance of a scalar 
electrostatic potential, which has to be added to the ef- 
fective Hamiltonian ([7]). 
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III. SCATTERING APPROACH 

In this Section we formulate the scattering approach 
to transport through a deformed graphene sample with 
metallic leads. We take advantage of the effective single- 
particle Hamiltonian 



H = -iaV 



V, 



(10) 



where the fictitious vector potential Ag is related to the 
strain tensor, u{R}, by means of the relation ([9]), and 
the scalar field V'(-R) describes strain-induced and exter- 
nal electrostatic potentials. In most of the intermediate 
expressions we let hv = 1 for simplicity. 

The metal leads are modeled by letting V = — C/ioad 
for a; < and x > L,^^ in Eq. (ITUl) . The vector potential 
is assumed to be zero in the leads. The width of the 
sample in y direction is denoted as W . The scattering 
off the metal leads are fully taken into account in the 
subsequent analysis. 

The rectangular sample geometry makes it convenient 
to employ the Fourier transform in y, 



R 



q 



(11) 



where q is the quasiparticle momentum component par- 
allel to the graphene-metal interface. The momentum 
takes on the quantized values, q = which depend 
on the boundary conditions in y direction, for example, 
qn ~ 2Tm/W with integer n for periodic boundary con- 
ditions. In the limit W ^ L, which we assume below, 
the particular type of the boundary conditions is not im- 
portant. 

We also restrict our consideration to small energies, 
\e\ Ulead, where the energy e is measured with respect 
to the Dirac point. In this approximation we derive the 
equation on the transfer matrix in the form— 



dM 

dx 



a^{q + T^Aey) +i<7z{V-e) +iT^Ag^ M, (12) 



where we introduced matrix notation in Fourier (channel) 
space, e.g. 



V„,. 



w 



dye'('J--''-^yV{x,y), (13) 



and qnm — Snmqn- The transfer matrix has the following 
structure in cr-space 



M = 



1/t't 



l/i 



(14) 



where f(f') for x = L is the matrix of reflection am- 
plitudes for quasiparticles entering the sample from the 
left(right) lead. The matrices t and t' contain the corre- 
sponding transmission amplitudes. Then, the Landauer 
formula for the conductance can be cast in the following 
form 



G 



2e2 



-TrttT 



2e2 



-Tr 



(15) 



where the trace is referred to the valley and channel space 
and the symbol TWii stands for the 11 block of the trans- 
fer matrix in cr-space. 

The Hamiltonian (|10|) with vanishing electrostatic po- 
tential, V = Q, obeys the chiral symmetry a^Hoz = —H, 
which is responsible for a non-Abelian Aharonov-Casher 
gauge invariance at zero energy. 



37 



(16) 



The spatially dependent phases and x can be cho- 
sen in such a way that the zero-energy spectral equation 
i?^o = is reduced to the Dirac equation, — icrV^'o — 0, 
with zero vector potential. This gauge transformation 
can be applied in the scattering approach in order to 
demonstrate that the presence of arbitrary vector po- 
tential has no effect on charge transport at the Dirac 
point as far as the contribution from edge states can be 
disregarded . '^^'•^^ 



IV. ID DIRAC-KRONIG-PENNEY MODEL 
A. Transport 

One-dimensional modulations of strain were realized 
experimentally in suspended graphene films using the 
remarkably large and negative thermal expansion of 
graphenCfiS Motivated by these experiments we calculate 
the transport properties of the one-dimensional Dirac- 
Kronig-Penney model, which has been introduced earlier 
by several authors.—^— Below we consider a general form 
of the model, where the variation of both electrostatic as 
well as vector potentials is included. Using the scattering 
approach formulated in the Section lllll we provide simple 
analytical solutions for transport and density of states in 
this model, which complement previous theoretical stud- 
ies of Dirac Fermions in periodic potentials.—"— 

A graphene sample with one-dimensional modulation 
of strain is characterized only by the component Mxa;(ic) 
of the strain tensor, which depends solely on the x- 
coordinate. From Eq. ([9]) , the components of the pseudo- 
vector potential take the form 



Aex = cos 36', Agy = 



■ sin 36* . 



(17) 



In addition, the strain induces a spatial variation of elec- 
trostatic potential, V = V{x). Its relation to the strain 
tensor is, however, complicated by screening effects. 

Since the potentials depend only on x, the transverse 
momentum q is conserved. In this case the matrices V, 
Agx and Agy in Eq. p2)) are diagonal in channel space. 
Since the considered potentials also do not couple the 
valleys, the scattering problem is reduced to the solution 
of 2 X 2 matrix equation. 



dM 

dx 



[cTx{q + Agy)-ia,ie-V) + iAgx]M. (18) 
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The a;-coniponent of the vector potential, Ag^, enters the 
equation in a trivial way and can be excluded by the 
gauge transformation M Aiexp{i Agj:dx'), which 
does not affect any observable. The first and the most 
trivial consequence of this transformation is that the one- 
dimensional strain modulations in x direction have no 
effect on transport in the zigzag direction (6 — 27m/ 3 
with integer n) provided the effect of V is ignored. 

For other angles the one-dimensional strain modula- 
tions lead to the appearance of new minima in the gate 
voltage dependence of the conductivity of the graphene 
sample. The numerical solution of Eq. (IT^ suggests 
that a periodic vector potential induces much more pro- 
nounced minima than a periodic scalar potential of equiv- 
alent amplitude. The weaker effect of V can be associated 
with Klein tunnelling, which leads to the suppression of 
pseudo-gaps in the spectrum of the superstructure. 

In order to describe the effect of periodic potentials 
analytically we introduce the one-dimensional Dirac- 
Kronig-Penney model, which is characterized by the vec- 
tor potential A given by Eq. ^ with 

uU^) = -uo/2 + uo^e(2;-2^j)e((2j+lK-a;), (19) 
j 

where the scale £ stands for half the period of the super- 
lattice, Q{x) is the Heaviside step function, and the di- 
mensionless parameter uq specifies the amplitude of the 
strain modulation. The vector potential introduced by 
Eqs. (|17ll9p corresponds to a strain field which is smooth 
on atomic scale but changes abruptly on distances smaller 
than the Dirac quasiparticle wave length, Tiv/e. The pe- 
riodic electrostatic potential is introduced in a similar 
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FIG. 2: (Color online) Conductance and Fano factor calcu- 
lated from Eq. (|28p for finite systems comprising A'' = 10 
and A'' = 50 supercells are compared with the averaged ex- 
pressions to be derived in Eqs. (|33l34p . which correspond to 
the limit N ^ oo and ignore contributions from evanescent 
modes. The plots are for graphene with strain modulations, 
rgx — 1. In a small vicinity of the Dirac point, e <^ hv/L, 
(shown at the insets) the transport is insensitive to strain due 
to the extended gauge invariance (|16|) and is dominated by the 
evanescent modes. 



V{x) = -Vo/2 + VoJ2 - 2^J) Q{{2j + l)i-x), (20) 
j 

where Vq is the amplitude of the modulation. 

For the sake of simplicity we consider a sample of 
length L = 2£N with N being an integer. In this case 
the solution to Eq. (fT8|) satisfies M{L) — Mq , where 
A^o is the transfer matrix corresponding to the size of 
the supercell, 2£. Disregarding the a;-component of the 
vector potential, Ag^, which has been argued above to 
have no effect on transport, we write 



Mo = M+M- 



(21) 



where we take advantage of the definitions 



q±=q± {Puo/2a) sin 36*, e± = e ± Vb/2. (22) 

The eigenvalues of the transfer matrix, Mq, are conve- 
niently parameterized as 



where we introduced the real function 

q+q--e+e 



A = cosfc+£ cosA;_£-|-- 



(23) 

s\nk+l sink-l, (24) 



with k± = (£j_— g^)^/^ the x-component of the momenta. 
The wave number fco becomes imaginary for some values 
of q and e indicating the appearance of pseudo-gaps in the 
superstructure. In the considered Dirac-Kronig-Penney 
model the value of fco is nothing but the a;-component of 
the quasi-momentum in the reciprocal space associated 
with the superstructure. 

In order to calculate the element Mu of the full 
transfer-matrix A4{L) we use the Chebyshev identity to 
calculate the A'^-th power of an unimodular matri^s^ 

A^jv A^osinfcoL- sinfco(L-2£) 
sm 2fco " 

The matrix element (A^o)ii = A — i?7 is readily deter- 
mined from Eqs. (|21l25p with 

?; = - — sin k-£ cos k+£ + sin k+i cos k-£. (26) 

fc_ fc-|_ 

Using the fact that both A and rj are real functions of 
the energy e, and the conserved momentum component 
q, we obtain from Eqs. (|14l25p the exact transmission 
probabilities for each value of q, 



cos knL 



1- A2 



sin fcgL 



(27) 
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FIG. 3: The density plot for the averaged transmission spec- 



tra, Tq{e), calculated from Eq. (|32|) for the Dirac-Kronig- 
Penney model with strain modulation rg^ = 0.1 (top left), 
= 2.0 (bottom left) and electrostatic field modulation 
ro = 0.1 (top right), ro = 2.0 (bottom right). 



This expression is reduced to the well-known result for a 
purely ballistic systemj^^ i.e. for uq — Vq — {), hy the 
substitution -q = (e/fc)sin2fc£ and A = cos2fc£. The 
transmission probabilities (|27|) determine the energy- 
dependent transport quantities 



G 



(28) 



where G is the Landauer conductance and F is the Fano 
factor for the shot noise. In the limit W ^ L the sum- 
mation over the scattering channels can be replaced by 
the integration, J2q ~^ (W/2t:) J dq. 

More generally one can define the cumulant generating 
function for the transport. 



^(x)-^ln(l-T, + e>^r, 



(29) 



The dimensionless cumulants c„ = Um^_^o (?"-^/9x" de- 
termine the so-called full counting statistics of the charge 
transport. The conductance and the Fano factor are 
given by G = (4e^//i)ci and F = c^/ci, respectively. 

The strength of strain-induced and electrostatic poten- 
tials in the Dirac-Kronig-Penney model (|19I20I) is char- 
acterized by the dimensionless parameters 

rg^ = (2uo£/a) sin 36*, ro = IVo/fiv, (30) 

respectively. In Fig. [2] the conductance and the Fano 
factor in graphene with periodic modulations of strain 
calculated from Eq. (|28l) are plotted for systems with 
finite lengths L = 20i and L = lOOi for re^ = 1. 
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FIG. 4: (Color online) The averaged conductance and Fano 
factor calculated from Eqs. p3l34p for the one-dimensional 
Dirac-Kronig-Penney model with different field strengths 
(|30|l . The effect of the scalar potential is suppressed due to 
Klein tunnelling through potential barriers. 



The energy dependence of conductance and Fano- 
factor reveal fast Fabry-Perot oscillations on the 
scale hv/L. From a physics point of view, these oscil- 
lations originate from multiple reflections of propagating 
modes (real fco) at the metal-graphene interfaces. 

The channels with imaginary values of ko (evanes- 
cent modes or metal-induced states) also contribute to 
transport. Even though the individual contribution 
of each evanescent mode is exponentially small, Tq ~ 
exp(— 2LImfco), their combined effect becomes essential 
for energies in a vicinity of band-edges. The role of 
evanescent modes is especially important at the Dirac 
point due to the absence of propagating modes. Indeed, 
the conductance and the Fano-factor determined by (pS)) 
take on the values G = Ae^W/nhL and F = 1/3 for 
e <C hv/L irrespective of the vector potential. This uni- 
versality is due to the extended gauge invariance (fT6| . 

For = L/2£ > 1, the amplitude of the Fabry-Perot 
oscillations in G and F decreases (provided W ^ L) and 
the relative contribution of evanescent modes becomes 
less important. Taking the limit A'^ — ^ c» is equivalent to 
ignoring the imaginary values of fco and averaging over 
the rapid phase, koL. This approximation has been used 
e.g. in Ref. [s^ to obtain the full counting statistics of 
few-layer graphene. 

The averaged generating function takes the form 



^(X) = 2^1n 



eX - 1 



(31) 
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FIG. 5: (Color online) The averaged conductance of Eq. (|33p 
(solid) is compared with an estimate (|35[) based on the mode- 
counting argument (dashed) for the Dirac-Kronig-Penney 
model with re^ = 0.3, ro = 0. The upper panel shows the 
corresponding transmission spectrum. 

where we introduce the mean transmission probability 
T, - |r;|-iReVl-A2. (32) 

The averaged transport quantities do not depend on the 
system size L and reveal no Fabry-Perot oscillations. For 
a ballistic sample, uq = Vb = 0, one finds — k/\e\, 

where k — ~ is the a;-component of the momen- 
tum. 

From the generating function pip we readily find the 
averaged conductance and noise. The conductance is 
given by the Landauer formula, 

G=— 5]T„ (33) 

while the Fano-factor is related to the averaged transmis- 
sion probabilities in a less evident way 

F^— ^. 34) 

The exact and averaged transmission probabilities 
(|27I32I) together with the corresponding expressions for 
the full counting statistics (|29l34p provide the complete 
analytical description of the transport properties in the 
ID-Dirac-Kronig-Penney model with scalar and vector 
potentials of arbitrary strength. 
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FIG. 6: (Color online) The averaged density of states obtained 
from Eq. (|43|) for the Dirac-Kronig-Penney model with differ- 
ent strengths of the periodic potentials. 

The dependence of the transmission coefficient on the 
transversal momentum component, g, is called the trans- 
mission spectrum. In Fig. [3] we plot the transmission 
spectra obtained from Eq. ((32|) for different strengths of 
the scalar and vector potentials. 

The transmission spectra shown in Fig. [3] are qualita- 
tively different for weak, r^x ^ 1, and strong, tq^ 3> 1, 
strain modulations. It has to be stressed that the latter 
regime is well within the experimentally accessible range 
of parameters. Indeed, for an achievable strain modula- 
tion of 0.2% (uo = 0.002) and the period i — 70 nm one 
finds re^ « 2. 

Neither the modulated strain nor the electrostatic field 
does open up a full band gap for any values of the pa- 
rameters or rg. Nevertheless, very large pseudo-gaps 
are generated by strain modulations with tq^^ ^ 1. One 
can see from the lower left panel in Fig. |3] that the charge 
transport in this regime is suppressed for almost all di- 
rections of the momenta (different q values) in a wide 
energy range. In contrast, the pseudo-gaps do not open 
at normal incidence (q — 0) in the Dirac-Kronig-Penney 
model with periodic electrostatic potential due to the 
Klein-tunnelling of quasiparticles through electrostatic 
barriers*^ 

The effects of strain and electrostatic field modulations 
on the conductance and noise are compared in Fig. 2] 
The most prominent feature of these plots are the dips in 
conductivity centered around the energies £„ = wnhvlll 
with n = 1,2,3,.... The physical mechanism, which 
is responsible for the appearance of the pseudo-gaps, is 
equivalent to that of the parametric resonancCfS^ 

For 3> 1, the conductance is suppressed in wider 
energy intervals around e = e„. For weak strain or po- 
tential amplitude, tq^ ^ 1 or ^ 1, significant pseudo- 
gaps are only located around the energy values e„ with 
odd values of n, since only odd Fourier components of the 
potentials (I19l20p exist. For a weak harmonic potential, 
a significant pseudo-gap would only arise around e — £\. 

The areas of vanishing transmission probability in 
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Fig. [3] (top left) can be understood on the basis of pertur- 
bation theory as presented in Appendix |X] for smah am- 
pHtudes of periodic potentials. The pseudo-gap emerging 
due to a weak periodic strain is described by the func- 
tions £g-{q) and eg-^-{q) given in the first row of Tab. ID 
These functions are shown with dashed lines in the upper 
panel of Fig. [S] 

Inside the pseudo-gap, i.e. for £g-{q) < £ < eg+{q), 
the transmission coefficient is exponentially small, while 



it is close to the ballistic value, T 



■ball 



outside the pseudo-gap (one can disregard the change in 
the resistance of the metal-graphene interface due to the 
weak periodic potentials). Note further, that the pseudo- 
gap is located near q — for energies near the dip such 
that « 1 in this region. We can, therefore, estimate 
the conductance in the vicinity of the lowest dip (e = Si) 
by subtracting the contribution of ballistically propagat- 
ing modes inside the pseudo-gap from the ballistic result 



G™c = Go - ^ ^ eie-Sg-iq)) e{eg+{q)-e), (35) 



where Q{x) is the Heaviside step function, Gq — e^We/h, 
and the limit W\e\ » 1 is assumed. The generalisation of 
Eq. (j35[) around the higher resonant energies, e = e„, n = 
3,5,..., is straightforward. 

It is shown in Fig. [5] for rg^ = 0.3 that the re- 
sult of Eq. ([35|) agrees with the averaged conductance 
calculated from Eq. ([55)1 . The band-structure analy- 
sis (|35IA5p predicts characteristic dips in the conduc- 
tance at £ = mrhv/2£,n = 1,3,5... of the depth 
SG = {8e'^W/h7T)y^\An\mT/2£, where A„ stands for n- 
th Fourier-component of the vector potential. For the 
Dirac-Kronig-Penney model one finds A„ oc 1/n hence 
the value of the conductance dip does not depend on n. 
We note that the validity of Eq. ([35)) is restricted to weak 
potentials. For rg^^VQ > 1, the pseudo-gaps overlap and 
the perturbation theory of Appendix |A] is no longer ap- 
plicable. 



Then, the integrated density of states per unit volume is 
given by 



q 



(37) 



where the factor 4 takes into account the spin and val- 
ley degeneracy. 

The transmission amplitude is readily obtained from 
Eqs. (ITinSi) as 



^ _ 1 
' cos koL + sin koL ' 



sin2fco^' 



(38) 



where the wave-number kg and the real quantity rj are 
determined from the expressions ([23]) and (|26p . respec- 
tively. Therefore, the partial density of states can be 
calculated exactly as 



L(q{dko/de) + sinkoL coskoL {d(q/de) 
' TT (cos^ koL + Q sin^ fco^) 



(39) 



Unlike the conductance or the shot noise ([28| . the den- 
sity of states depends on the value of C/iead, because the 
spectrum in a vicinity of the band-edges is dominated by 
the metal- induced states (evanescent modes). 

The metal proximity effect^^ can be seen already for 
purely ballistic system near the Dirac point. Indeed, for 
uq = Vq = 0, one finds 



hall 



(e) 



Le^ — {q^ /k) cos kL sin kL 
7r(e^ — q^ cos^ kL) ' 



(40) 



where k = \/ — q^- At zero energy the ballistic result 



H0| is reduced to ^^{0) = ta.nh{qL)/Trq. Therefore, 
the density of states at e = in Eq. (I57|) acquires a loga- 
rithmic divergency. This divergency is regularized by the 
largest available transversal momentum (?niax, which is 
simply equal to the Fermi-momentum in the metal lead, 
<Zmax — Uicad/^v- Thus, the density of states in a close 
vicinity of the Dirac point, e hv/L, is given by 



^'^""(O) 



■In 



UlcadL 



(41) 



B. Density of states 

Let us now extend the solution of the Dirac-Kronig- 
Penney model to the density of states. For purely ballis- 
tic system with metal leads, both the local and the inte- 
grated density has been found in Ref. using a Green's 
function approach. Below we take advantage of an alter- 
native route and relate the partial density of states in the 
channel q to the corresponding transmission amplitude, 
tq, by the well-known formulai^ 



in agreement with Ref. i56 . Similar logarithmic depen- 
dence of the density of states on f/ioad takes place in 
the Dirac-Kronig-Penney model for energies inside the 
pseudo-gaps. 

In full analogy with the averaged generating function 
([5T|) we can introduce the averaged density of states, p, 
which corresponds to the limit iy — > oo. In this limit we 
disregard the contribution of the metal-induced states by 
projecting on the real values of the momentum fcg and 
average the result of Eq. (l40l) over the rapid phase /cqL. 
This procedure leads to the simple result 



Vq ~ Im 

TT 



d In tn 



de 



(36) 



L 

TT 



dRc fcn 



ds 



(42) 
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FIG. 7: (Color online) The contour plot of the dispersion re- 
lation for the Dirac-Kronig-Penney model (|19l20p with strain 
modulation rex = 0.3 (top left), re^ = 2.0 (bottom left) and 
electrostatic field modulation ro = 0.3 (top right), ro = 2.0 
(bottom right) at energy e — 3.3hv/2£. 



hence, in the limit W ^ L, the mean density is given by 



P^-K dq 



dRe fco 



de 



(43) 



In the ballistic limit, uq — Vq — 0, Eq. (|43p yields the den- 
sity of states of the clean graphene p*^^'' — 2\e\ / nh'^ v"^ . In 
Fig. |6] we plot the averaged density of states calculated 
from Eq. (|43l) for different strengths of the periodic po- 
tentials. 



C. Beyond the Dirac-Kronig-Penney model 

The results of the previous subsections are readily gen- 
eralized to a model with arbitrary periodic variation of 
strain and electrostatic potential in x direction. The 
analysis of the model is reduced to the calculation of 
the transfer matrix, A4o, which corresponds to the wave 
propagation over the distance 2£, the period of the po- 
tential. In this generalized model, both the exact and the 
averaged full counting statistics as well as the density of 
states are still given by the expressions (|29I31I39I42[) with 



(44) 



and the a;-component, fcoj of the quasi- momentum is 
related to A by Eq. 1^^. Thus, for a periodic potential 
of a general type, the full solution of the problem is re- 
duced to the straightforward numerical evaluation of the 
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ky ky 



FIG. 8: (Color online) Typical band gap boundaries. In all 
plots n, b was set to 1. 



functions A and rj. Note, that the exact analytical expres- 
sions (j24l26l) are restricted to the Dirac-Kronig-Penney 
model and do not apply generally. 

We have used Eqs. (|44|) to calculate the energy- 
dependent conductance for different amplitudes of the 
periodic strain in the harmonic potential, u^x = 
1*0 sin(7ra;/.^). In this case, the higher pseudo-gaps are 
found to be suppressed for weak potentials re^ , tq ^ 1 as 
compared to the Dirac-Kronig-Penney model. The low- 
est conductance minima around e = is essentially the 
same in both models. The models become even more sim- 
ilar with increasing potential strength. For rg^ ,rQ }^ 1 
the transmission spectra of the single-harmonic model, 
becomes almost equivalent to those of the Dirac-Kronig- 
Penney model for all energies. 



D. Transport in y direction 

We have demonstrated that periodic and x-dependent 
scalar and vector potentials modify transport in x- 
direction due to the appearance of pseudo-gaps. Let us 
now argue that such potentials hardly affect transport 
in y direction, i.e. in the direction parallel to the cqui- 
potential lines. 

In Fig. [7| we plot the Fermi surface slightly above the 
energy value e = Eq using the exact dispersion relation, 
e = e{ko,q), obtained from Eqs. (|23l24p for the Dirac- 
Kronig-Penney model. We have seen that the transport 
in X direction is dominated by the modes with the mo- 
menta q € (— For some values of q in this interval 
one cannot find any propagating states, which correspond 
to the real values of fcg. This can be seen as the forma- 
tion of pseudo-gaps in the transmission spectra shown in 
Fig. [21 which is especially strong for r^a;, ro > 1. In con- 
trast, for each real value of fco one always finds real values 
of q. Therefore, the transport in j/-direction is not af- 
fected by the formation of the pseudo-gaps, and depends 
only on the details of scattering at the metal-graphene 
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interfaces. The latter effect is, liowever, weak and falls 
beyond the scope of our consideration. We therefore con- 
clude that one-dimensional superlattices have a major 
impact only on the transport properties along the direc- 
tion of their periodicity. 



V. CONCLUSION 

The Dirac-Kronig-Penney model has been applied in 
order to study phase-coherent charge transport in peri- 
odically strained graphene samples. Using the exact re- 
lations for the band-structure in the considered superlat- 
tice, we calculate the density of states and the full count- 
ing statistics for charge transport. The exact quantities 
are simplified in the thermodynamic limit by neglecting 
the contribution of evanescent modes and averaging over 
the Fabry-Perot oscillations. The conductance is found 
to be suppressed at energies corresponding to the posi- 
tions of the pseudo-gaps in the modified spectra of the 
superlattice. The periodic strain and electrostatic po- 
tentials are argued to have the largest effect on transport 
in the direction perpendicular to the equipotential lines. 
The vector potentials are shown to play a greater role in 
confining the Dirac quasi-particles due to the suppression 
of Klein tunnelling. 
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Appendix A: Band structure for Dirac particles in a 
superlattice 

In this appendix we perturbatively calculate the band 
structure of a graphene sample described by the single- 
valley Dirac equation 



where 
V(r) 



[-io- ■ V + V(r)] *(r) = e*(r), 



Vir) + VAr) Ag.,ir) - iAeyir) 
Ag.,{r) + lAgyir) V{r)~V,{r) 



(Al) 



(A2) 



is a general Id-periodic potential that is smooth on the 
scale of the graphene lattice constant a. To do so, we gen- 
eralize the standard methods of elementary band struc- 



ture theory to the Dirac case with its non-trivial ma- 
trix structure. We find that for purely scalar potentials, 
this additional structure prevents the creation of a gap 
at ky — thus making Ag and Vz more relevant for 
the transport properties; this fact can be related to the 
chirality of the quasi-particles in graphene and the phe- 
nomenon of Klein tunnelling. 

Let b — [bx,by)^ span the Id-super-lattice (so that 
V{r + nb) ~V{r)) and g be the corresponding reciprocal 
vector fixed by the conditions g ■ b = 2tt and g \\ b. 
Making the Bloch ansatz 



*(r) 



E 



il){k — ng) e 



i{h~-ng)r 



(A3) 



for the wave function and writing the periodic potential 
as a Fourier series, V(r) = X^n "^'n e"'^'', we can write 
down the central equation 



— £ k* — Tig 
k — ng — £ 



i>{k-ng)+^Vjn^n-4>{k''mg) ^ 0, 

in 

(A4) 

where 'k = kx + iky and g — gx + igy The central equa- 
tion (|A4[) represents an infinite set of linear equations for 
the coefficients V'l which can in general only be solved 
numerically. The solution, however, simplifies for suffi- 
ciently weak periodic potentials. 

For k close to a degeneracy, \k\ — \k — ng\ ^ e, we can 
neglect all but the two coefficients tp{k) and 4'{f^ ~ "-9) 
in Eq. (|A4p . The resulting linear system has the simple 
form 



Ag 



■14 

iAg 



A 



V* -f 

An„ + iA^ 



Al 



By 



V* 



-v:' 



9x - iAg 

V-Vz 
k* 



kn 



(A5) 

Here, we introduced the abbreviation k„ = k — ng and 
omitted the Fourier index n of the potentials for brevity. 
The solution of Eq. (|A5|) is particularly compact when 
only one of the coefficients Agxn, Agyn,Vn or Vzn is dif- 
ferent from zero. In this case we obtain 



|k„ 



± iV(|kP-|k„P)2+4|C/„|2|2p-ng|2, (A6) 

where the pair of quantities (C/„,p) take on one of the fol- 
lowing values {Axn.kx), {Ayn,-iky), (U„,fc), or (14„,0). 
We compare the energy values of the forbidden zone's 
upper/lower boundary, £g±, as well as the width of the 
gap in forward direction for the different types of periodic 
potentials in Tab. HI To calculate £g± we have chosen k 
in Eq. (jA6[) such, that k-b/b = n|g|/2. Typical band gap 
boundaries are shown in Fig. [S] 
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A{ky = 0) 


Aexn 7^ 


^k:i + \^f±\Ae.„\ 





Aeyn ¥= 




'2\Aey„\ 




^Jk^y + \^\' + \Vn\^±\ky\\Vn\ 







^/k^y + m' + \V.n\'±2mV.n\ 


2\V.n\ 



TABLE I: Energy value of the forbidden zone's upper/lower 
boundary, as well as the width of the gap in forward direction 
for the diflferent types of periodic potentials. In the first row, 
b was chosen parallel to By, for the remaining rows b \\ e^. 
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